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ABSTRACT 

Arecibo observations of the conal triple pulsar B1918+19 at 0.327- and 1.4-GHz 
are used to analyse its subpulse behaviour in detail. We confirm the presence of 
three distinct drift modes (A,B,C) plus a disordered mode (N) and show that they 
follow one another in specific cycles. Interpreting the pulsar's profile as resulting 
from a sightline traverse which cuts across an outer cone and tangentially grazes 
an inner cone, we demonstrate that the phase modulation of the inner cone is 
locked to the amplitude modulation of the outer cone in all the drift modes. The 
9% nulls are found to be largely confined to the dominant B and N modes, and, in 
the N mode, create alternating bunches of nulls and emission in a quasi-periodic 
manner with an averaged fluctuation rate of about 12 rotation periods (Pi). We 
explore the assumption that the apparent drift is the first alias of a faster drift 
of subbeams equally spaced around the cones. This is shown to imply that all 
modes A, B and C have a common circulation time of 12 Pi and differ only in 
the number of subbeams. This timescale is on the same order as predicted by the 
classic E X B drift and also coincides with the N-mode modulation. We therefore 
arrive at a picture where the circulation speed remains roughly invariant while 
the subbeams progressively diminish in number from modes A to B to C, and 
are then re-established during the N mode. We suggest that aliasing combined 
with subbeam loss may be responsible for apparently dramatic changes in drift 
rates in other pulsars. 
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1 INTRODUCTION 

Pulsar B1918+19 is a typical 'slow' radio pulsar with a 
stellar-rotation period (Pi) of 0.82 s and a spin-down age 
of 15 Myr. It was one of three apparently conal pulsars 
with three-peaked profiles whose pulse-sequence (here- 
after PS) behaviour was studied at 430 MHz by Hank- 
ins & Wolszczan (1987; henceforth HW). These authors 
identified four distinct modes of emission labelled A, B, 
C, and N — a nomenclature that will also be used here 
for continuity. The first three modes exhibited drifting 
subpulses which, together with the fourth "disordered" 
mode and embedded null sequences, were found to inter- 
act in specific sequences. 
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In this paper we exploit the improved signal-to-noise 
ratio (S/N) of our observations to analyse this pulsar's 
drift modulation and nulls more thoroughly. A full sta- 
tistical assessment is given for each mode and type of 
modal sequence. Recent studies of other similar pulsars 
have found that nulling is sometimes confined to a specific 
mode of emission (Redman ei al. 2005) or found to have 
a periodic nature (Rankin & Wright 2007; Herfindal & 
Rankin 2007, 2009) and we establish here to what degree 
this is true of B1918-f 19. 

To place our results in a physical, or at least geo- 
metric, framework, we consider the context of a rotating 
subbeam-carousel system (RS). This has often proved 
successful in describing drifting subpulses in conal pul- 
sars {e.g., Deshpande & Rankin 1999, 2001; hereafter DR; 
Battacharyya et al. 2007), and it is now interesting to ex- 
plore whether the more complex multiple discrete mode 
patterns of B1918-I-19 together with their 'sequence rules' 
can be made compatible with such a model. 
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Table 1. Band, resolution, and length of Arecibo observations. 



Band 


MJD 




Resolution 


Length 


(GHz) 




(d 


agrees/sample) 


(pulses) 


1.1-1.7 


52735 




0.22 


1023 


0.327 


52942 




0.22 


1395 


0.327 


53778 




0.39 


3946 


1.1-1.7 


54541 




0.36 


731 



The first t\vo observations \vcrc seriously eorrupted by interference. 

In §2 we describe our observations, and in §3 we 
present what can be inferred about the pulsar's basic 
emission geometry from its profile and polarization. §4 
introduces the multiple features in its fluctuation spec- 
tra, and §5 discusses the profiles and phase patterns of 
the several modes. §6 gives an analysis of nulling statis- 
tics, and §7 discusses mode sequencing and transitions. 
In §8 we discuss the possibility of aliasing, and §9 pro- 
vides an interpretation of the pulsar's three drift modes 
in terms of the rotating carousel-beam model, and this 
is checked against the observationally-derived geometry 
of the pulsar in §10. §11 summarises the results and §12 
draws conclusions. 



2 OBSERVATIONS 

The observations were carried out using the 305-meter 
Arecibo Telescope in Puerto Rico. All of the obser- 
vations used the upgraded instrument with its Grego- 
rian feed system, 327-MHz (P band) or 1400-MHz (L 
band) receivers, and Wideband Arecibo Pulsar Proces- 
sors (WAPlQ). The ACFs and CCFs of the channel volt- 
ages produced by receivers connected to orthogonal lin- 
early (circularly, after 2004 October 11 at P band) polar- 
ized feeds were 3-level sampled. Upon Fourier transform- 
ing, some 64 or more channels were synthesized across 
a 25-MHz bandpass with about a milliperiod sampling 
time. Each of the Stokes parameters was corrected for 
interstellar Faraday rotation, various instrumental polar- 
ization effects, and dispersion. The date, resolution, and 
the length of the observations are listed in Table [2] 



3 PROFILE GEOMETRY 

The B1918+19 polarization profile and PPA histogram, 
shown in Figure [l] consists essentially of four compo- 
nents, a distinct leading component, a central "main" 
component exhibiting a clear shoulder, plus a weak trail- 
ing component some 35° of longitude after the peak. Its 
profile has therefore been placed in the relatively un- 
common conal quadruple (cQ) category (Rankin 1993; 
ET VI), and we now know that some other such stars 
also have highly asymmetric profiles {e.g., J1819-f 1305; 
see Rankin & Wright 2008). The available profiles of 
B1918+19 {i.e., Gould & Lyne 1998; Hankins & Rankin 
2008) suggest that it becomes more symmetric at higher 
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Figure 1. Average profile of B1918-I-19 computed from the 
long 327-MHz observation. Notice the leading and compound 
primary components and the relatively weak trailing compo- 
nent. The top panel gives the the total intensity (Stokes /; solid 
curve), the total linear (Stokes L {=\/ Q^ + U^)\ dashed), 
and the circular polarization (Stokes V; dotted). The lower 
panel PPA (=^ tan~^((7/Q)) histogram corresponds to those 
samples having PPA errors smaller than 10° with the average 
PPA traverse overplotted. The longitude origin is taken about 
halfway between the outer conal component pair. 



frequencies, but no published profile at a frequency higher 
than 1.6-GHz seems to exist. The relative intensity of the 
first component at meter wavelengths also seems to vary 
considerably, but as we will see below, this is probably 
due to the specific population of modes present in the 
average. 

The aggregate linear polarization L never exceeds 
about 10% in Fig. [1] but the geometric polarization an- 
gle (hereafter PPA) sweep rate can be reliably traced 
once the effects of two modal "flips" are resolved. The 
PPA then makes a total traverse of about -120°, and the 
central rate R is some -3.2°/°. Despite the then limited 
information, B1918-|-19's basic emission geometry seems 
to have been computed almost correctly in ET VI: its 
respective inner and outer cones have half-power, 1-GHz 
widths of 20.5 and 51.3°, respectively. Its magnetic lat- 
itude a and sightline impact angle /? are then 14° and 
-4.3° for a poleward traverse across its emission cones (or 
11° and -1-3.4° for an equatorward traverse, although this 
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Figure 2. Two single pulse plots, each of 340 consecutive 
pulses from the same observation of B1918+19, separated by 
about 1000 pulses. Left: several sequences of mode A trans- 
forming into mode B interposed by the null-rich irregular N 
mode. Note the nulls creating an approximately 85-pulse pe- 
riodicity. Right; an extended sequence of the fast mode C in- 
terrupted by an extended null. Longitude scale in degrees. 
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model fits the inner/outer conal dimensions less success- 
fully). In either case, the sightline encounters the inner 
cone tangentially (with /3/p ^0.9, where p is the outer 
half-power conal radius) and the outer cone much more 
centrally (-0.65/0.52) — confirming the earlier qualitative 
interpretation of the star's double-conal emission geome- 
try. Moreover, the prominent polarization-modal activity 
on the profile edges seems to confirm that we are ob- 
serving emission from the full widths of the cones. Note, 
though, that the PPA distribution in the leading half 
of the central component consists of a set of "patches" 
rather than two PPA tracks. 



4 PULSE-SEQUENCE MODULATION 

The two extracts from our longer 327-MHz observa- 
tion presented in Figure [2] show the essential features 
of B1918-|-19's single pulse behaviour. Using the nomen- 
clature of HW, the 340 pulses on the left show several 
examples of a relatively slow drift pattern (A) changing 
rapidly into a faster mode (B), which is then followed by 
an irregular null-rich sequence (N). The same number of 
pulses in the right panel show a typical long sequence of 
the fastest mode (C) interrupted by a relatively long null. 
The exceptional S/N of the whole observation enables us 
to analyse these basic features in some detail. 

The longitude- and harmonic-resolved fluctuation 
spectra (hereafter Irf and hrf ) of the entire 327-MHz ob- 
servation are given in Figures [3] and [H In the Irf spectra 
four clear features are visible, in addition to a prominent 
low-frequency feature. The periodicities corresponding to 
the A, B, and C modes, which represent the vertical spac- 
ing of the driftbands and are conventionally denoted as 
P3, are very evident, with centroid values only slightly 
different from those of HW. The A-mode has a P3 of 
6.06 Pi (=1/0.165 c/Pi; 5.9 Pi in HW); the B-mode a 
P3 of 3.80 Pi (=1/0.263 c/Pi; 3.85 Pi in HW); and the 
C-mode a P3 of 2.45 Pi (=1/0.408 c/Pi; 2.5 Pi in HW). 
Substantial scatter about these values is visible in Fig. [3] 
(especially in the B-mode), but we will show below that 
they are often nearly coherent in particular sections of 
the PS. Note also that the periodicities corresponding 
to these three modes are found not only in the central 
component of the integrated profile, but also in the outer 
component pair. Significantly, the A and B modes modu- 
late both the leading and trailing components, although 
on the trailing edges at differing longitudes. The C mode 
modulation is also prominent in the leading component 
but invisible in the trailing one, where the intensity of 
this mode is anyway very weak|j 

The A-C features also appear prominently in the hrf 
spectrum of Fig. |4l but here almost exclusively at fluctu- 
ation frequencies between about 0.6 and 0.85 c/Pi. The 
main panel shows clearly that these periodicities repre- 
sent phase modulations corresponding to sets of drift- 
ing subpulses separated by P2 ~20°(=360°/18) and their 



2 Weltevrede ei al (2006, 2007) were unable to find modula- 
tion in this pulsar at 21 cms; whereas, at 92 cm they detected 
two drift features. See their fig. A12. 
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Figure 3. Longitude-resolved fluctuation-power spectra for 
the 327-MHz observation of B1918-I-19 on MJD 53967 (main 
panel, per the colour scale at the right), plotted against the 
total power profile (lefthand panel) with the integral spectrum 
in the bottom panel. Five different features are seen in these 
spectra: The three highest frequency features correspond to 
the A, B and C modes identified by HW with P3 values of 
5.8, 3.9 and 2.5 cycles/Pi. In addition two lower frequency 
periodicities are present, a well defined one corresponding to 
a P3 value of some 85 P\ and a broad response at about 12 Pi. 
Note that the 2. 5- Pi, 3.9-Pi and 5. 8- Pi features all modulate 
outer regions of the outer cone to some degree although the 
2.5-Pi feature is not detected in the trailing component. Fi- 
nally, the low frequency feature modulates the entire profile, 
whereas the 12-Pi modulation affects the central inner-cone 
region only. 



presence at frequencies over 0.5 indicate that the appar- 
ent drift is toward increasing longitudes. In itself the hrf 
spectrum cannot determine the true subpulse drift, i.e 
whether or not the observed drift pattern is aliased (in 
which case the intrinsic drift would be toward decreasing 
longitudes). This point is further discussed below. 

This new periodicity may well arise within the irreg- 
ular null-rich N-mode, which is known to have a peak in- 
tensity at a slightly later longitude than the other modes 
and to be characterised by a weak first component. This 
will be investigated in a later section, but it can be seen 
in the left-hand PS of Fig. [2]that intensity fluctuations on 
roughly this timescale occur between the A-B modal se- 
quences. It is interesting to note — not for the first time in 
pulsars with multiple drift patterns — that the periodici- 
ties of this pulsar appear to have "harmonic" properties: 
the new periodicity (12.0 Pi) is practically twice that of 
the A-mode (6.06 Pi), which in turn is about three-halves 
of the mean value of the B-mode (3.80 Pi) and five halves 
of the C-mode periodicity (2.45 Pi). 

Finally, the low-frequency feature in the two 
spectra — mentioned in passing by HW and attributed to 
"deep nulls" — has its peak before the profile peak and 
is prominent across nearly the entire profile. A striking 
example of how this feature is is generated in the PS can 
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Figure 4. Harmonic resolved fluctuation spectrum for the 
same observations as in Fig. |3] (main panel, per the colour 
scale at the right). Rather dramatically, this shows that HW's 
three features are all phase-modulated. The low frequency re- 
sponse, by contrast, is a mixture of phase and amplitude mod- 
ulation, and the 12-Pi response appears to represent mostly 
amplitude modulation. This distinction can also be seen in the 
main panel where the low frequency modulation is carried by 
smaller-numbered Fourier harmonics than HW's drift modula- 
tion, which seem to peak about number 18 or so. In the central 
panel the amplitude scale has been compressed by a factor of 
3 to accentuate the three phase-modulated periodicities. 



be seen in the left of Fig. [2l where the 340 pulses are 
divided into 4 roughly equal sections by the intrusion of 
the N mode. However, in the PS as a whole this rep- 
resents more of a characteristic timescale rather than a 
periodicity. In the hrf spectrum (Fig. |4]) the feature in- 
terestingly reveals itself to be an admixture of amplitude 
and phase modulation — the latter counter to that of the 
A-C modes and focussed on a narrow band of longitude 
preceding the profile peak by about 3°. (NB, a similar 
feature was found in J1819+1305 in an earlier paper.) 

As suggested by Fig. [2l these last two periodici- 
ties (12-Pi and 85-Pi) do appear to be null related. A 
number of the nulling pulsars exhibit such low frequency 
fluctuation-spectral features, and when the PS modula- 
tion is quelled by substituting a scaled-down profile for 
the pulses and zeroes for the nulls, the periodicities often 
persist (Herfindal & Rankin 2007, 2009). The latter study 
shows that this is so for B1918+19 as well; three different 
periodicities appear to be present in the pulsar's null se- 
quence: a seemingly harmonically-related pair at 85±14 
and 43±4 Pi as well as the roughly 12-Pi one. 



5 DISTINGUISHING THE MODES 

Using colour displays such as those in Fig. (2] we inspected 
the entire long 327-MHz PS to identify sections exhibit- 
ing the A, B and C modes, usually on the basis of their 
characteristic periodicity. Mode N was identified by the 
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Figure 5. Partial profiles for the several modes as well as 
the average profile: (top left) normalized to their relative in- 
tensities; (top right) leading component detail; and (bottom) 
scaled to the peak intensities. Total average profile (solid line) 
as well as the A (dashed), B (dash-dotted), C (dotted), and 
N (dash-triple dot line) modes. Notice that the leading com- 
ponent of the C mode is significantly larger than in any other 
mode, whereas its trailing component is much smaller. It also 
peaks about 1° earlier than the A mode, while the B mode 
peaks somewhere between. Similarly, the normalized profiles 
(bottom) clearly show the shift of the peak of the N mode as 
well as the curtailed trailing edge of the C mode. 
© 2008 RAS, MNRAS 000.[THT3l 
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absence of a clear drift, often accompanied by weakness 
in the first component. 

By far the most common mode was B (1896 pulses 
or 48%), followed by N (847 pulses or 21.5%), A (662 
pulses or 17%), and C (541 pulses or 14%). However, 
these statistics disguise the considerable discrepancies in 
the typical sequence lengths of the modes: while all C- 
mode pulses occurred in just 4 long null-free sequences, 
averaging 135 pulses in length, the 36 N-mode PSs aver- 
aged only 23.5 pulses, and the A- and B-mode PSs aver- 
aged 35 and 53 pulses, respectively. 

By segregating the pulses appropriately a partial- 
average profile was computed for each mode. Figure [5] 
(top left) shows the profiles weighted according to their 
relative average intensities. Note that the average inten- 
sities of the A- and C-mode pulses are clearly stronger 
than those of the B and N modes. This differs slightly 
from the conclusion of HW, who found, after selecting 
300 bright pulses from each mode (at 430 MHz), that the 
average intensities of all the modes were more compara- 
ble. Their result can be reconciled with ours if, as seems 
likely, the weaker pulses (including nulls) tend to occur 
in the B or N modes. Particularly striking in both ours 
and HW's results is the enhanced intensity of the leading 
component in the C mode, which is about double its av- 
erage strength in the other modes. The first component 
peak appears about 1° earlier in the C mode than in the 
A mode, with those of B and N modes falling in between 
(top right panel). 

By contrast. Fig. [S] (bottom) shows the same par- 
tial profiles, but now normalised to the amplitude of the 
total average profile peak. This exhibits the attenuated 
latter half of the C-mode profile: there is no "hump" on 
the trailing edge of its main component, and its weak 
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Figure 6. Drift-modal sequences folded at their respective 
P-i values: A (top left, pulses 2404-2458); B (top right, pulses 
2600-2800); and C (bottom, pulses 481-735). Note the strong 
similarity of all three patterns, including the relative phase po- 
sitions of the leading and trailing components. Thus, whether 
aliased or not, in each mode the pulsar is executing virtually 
the same pattern at a different speed. 



fourth component all but vanishes. The fourth compo- 
nent is strongest in the A mode and appears at a closer 
longitude to the main pulse than the fourth component 
of the B mode. The plot also enables comparison with a 
similar graph in HW (their fig. 4). The results are broadly 
compatible, but our profiles show a much smaller relative 
intensity in the region between the leading and main pro- 
file components, reflecting the much better resolution of 
our observations. 

HW also found that the peak in the N mode falls 
about 3° later than in other modes, whereas our results 
do indeed show a peak shift, but only of about 2° (Fig. [5] 
bottom left panel). A possible explanation for these dis- 



crepancies can be a shift in the centroid of the drifting 
pulses towards earlier longitudes as they become more 
intense. 

Although the modal partial profiles in Fig. [5] all have 
distinguishing properties, the underlying subpulse config- 
uration in the three drifting modes A, B and C is nearly 
identical. The three panels of Figure [6] display selected 
intervals folded at the appropriate P3 values for modes A 
(top left), B (top right) and C (bottom), respectively. In 
all three modes the drift bands appear to have a similar 
longitudinal separation (P2), a point also noted by HW. 
Most striking, however, is the phase-correlated behaviour 
of the leading and main-component emission patterns. 
The leading component does not drift but instead modu- 
lates its intensity at the same P3 as the drift in the main 
component (evident in the fluctuation spectra of Figs. [3] 
andg]). 

Note that the phase of the subpulse emission in the 
leading component is correlated with the drift patterns 
in the main components in a similar manner in all three 
modes. Although the emission level of the leading compo- 
nent varies from mode to mode, in each case it is strongest 
as or just before the drift band reappears at the leading 
edge of the main component. A similar correlation is also 
discernible between the main and trailing components. 
In the A mode (Fig. (H] top left) there is evidence of an 
emission peak at the longitude of the trailing component 
of the main pulse. Interestingly, this peak does not occur 
at the same phase as that in the leading component, but 
follows it by about 4.5 Pi or 0.65 in phase. There is also a 
hint of this peak at the same relative phase in the B-mode 
fold of Fig.[Sl top right, lending an inverse-reflection sym- 
metry to the modulation pattern [noted also between the 
two halves of the emission pattern of B0834-I-06 (Rankin 
& Wright 2007) and in B0818-41, where the phase shift 
is close to 0.5 ( [Bhattacharyya et al. 2009] )]. 

Thus B1918+19 joins the relatively short list of pul- 
sars both classifled as having double-cone structure and 
whose single-pulse patterns are phase-locked (B1237-I-25, 
B0826-34, B0818-41, B1039-19: see Srostlik & Rankin 
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Figure 7. Pulse-intensity histogram (solid curve) for the long 
327-MHz observation of B1918-I-19 and a corresponding off- 
pulse noise window (dashed curve). The vertical dotted line 
indicates a plausible null threshold at 0.18 </> (see text). 
Some 21 pulses were omitted from the distribution due to bad 
baselines. 



2005, Gupta et al. 2004, Bhattacharyya et al. 2007, 
2011 respectively). However, as pointed out by Bhat- 
tacharyya et al. (2009), there are as yet no examples 
where phase-locking is not present in double-cone pul- 
sars. In B1918-I-19 we can additionally see that the phase- 
locking survives multiple drift-mode changes. The simi- 
larity of the pattern in all three modes suggests that the 
emission processes are simply executing a very similar 
pattern at three different speeds and demonstrates that 
all components interact as a single physical system and 
points to a common geometry and a common emission 
process. 
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Figure 8. Histograms of null (top) and burst (bottom) 
lengths. 



6 QUASI-PERIODIC NULL PATTERNS 

The excellent S/N ratio of the PSs from the upgraded 
Arecibo instrument enabled us to study B1918-fl9's 
nulling properties in detail for the first time. Figure [7] 
gives the pulse-intensity histogram for the 327-MHz 3946- 
length PS. It shows a population of weak pulses around 
zero energy that is distinct from the low-energy tail of 
the primary intensity distribution. This suggested a con- 
servative definition of a null as having an intensity less 
than the minimum in the distribution, which occurred at 
0.18 </>. Using the indicated threshold at 0.18 </>, 
the null fraction was found to be ~8.9% or 350 pulses. 

A histogram of null lengths in Figure |8] (bottom) was 
computed to investigate their statistical properties. Its 
most obvious feature is the peak of 80 isolated null pulses. 
If the 350 nulls had appeared randomly among the 3946 
pulses with an 8.9% chance in each pulse, the expected 
number of isolated nulls would have been some 320. The 
non-randomness of the distribution is then obvious from 
the massive excess of null pulses occurring in bunches 
longer than one pulse. It is also reflected in the highly 
negative "over-clustered" Runs- Test value (see Redman 
& Rankin 2008). 

Visual study of the PSs (see Fig. [2]) strongly sug- 



gests that the nulls are closely associated with the N 
mode, occurring predominantly before, during, after, or 
within it — so it remains possible that nulls are appear- 
ing randomly, but confined to this mode of emission. If 
we postulate that all the nulls occur with this set of 847 
N-mode pulses — where its null fraction would then be 
41% (=350/847) — we could expect some 82 single nulls, 
very close indeed to the number that we observe. The 
null-length distribution is also statistically overweight in 
relatively long (>3) contiguous nulls. In particular, the 
observed 3 occurrences of 10 and 11 consecutive nulls 
remain highly improbable. It seems an inescapable con- 
clusion that the nulls in this pulsar do have an intrinsic 
bunching tendency such that once a null pulse occurs, the 
chance of a null in the next pulse is enhanced. 

We also noted that groups of null pulses often pre- 
cede and follow a typical 6-10 pulse burst of emission, 
even once during the B-mode in an apparent interruption 
of that mode. In Fig. [8] (top) a histogram of the emission- 
burst lengths indeed reveals a distinct peak around 6-10 
pulses, which can be attributed to this emission bunching. 
(Also visible is a statistical excess of 1-Pi bursts, reflect- 
ing emission "flickering" within sequences of nulls) . Thus 
nulls not only tend to occur in bunches, but those bunches 
appear in a quasi-periodic manner — obviously in the N 
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mode and often during transitions between modes (see 
next section). We have already noted (e.g., Fig. [5]) that 
the N-mode emission is characterised by weak leading- 
component emission and a slight shift in the main peak 
toward later longitudes (see Fig. [2|. These are precisely 
the properties shared by the unexpected 12-Pi feature 
in the fluctuation spectra of Figs. [3] and IH and we are 
therefore confident that the quasi-periodic nulls associ- 
ated with the N mode are responsible for that feature. 



slightly faster than the B mode, only to settle down to 
the typical B-mode drift rate within a few pulses. In the 
left panel of Fig. [2] it is also evident that each A — > B 
sequence has increasing then diminishing intensity, giv- 
ing the mode change the sense of evolution. This is not 
seen in the purely C-mode sequence shown in the right 
panel, suggesting that the two types of sequences are of 
different character, with the latter more stable than the 
first. 



7 MODAL SEQUENCES AND 
TRANSITIONS 

Studies of the few pulsars with multiple drift modes 
(Huguenin et al. 1972; Wright & Fowler 1981; Redman et 
al. 2005) have shown that the modes occur in sequences of 
increasing drift rate. HW found this to be true for pulsar 
B1918-f 19 also, and determined, more specifically, that it 
exhibits two characteristic modal sequence patterns: null 
^N->A^B^N^ null, and null -^ C ^ null. 

We found that these two patterns do indeed occur, 
and they are illustrated in Fig. [2] with several examples 
of the first type shown in the left-hand panel and the C- 
mode dominated pattern on the right. However, they are 
not the only patterns this pulsar exhibits. In particular, 
we noted frequent examples of the B and N modes being 
both preceded and followed by null bunches (see 33 , and 
there was one possible instance of a null — >■ B — >■ C — >■ 
null sequence. The A mode never seems to null and is 
invariably followed by B mode. It seems that in all cases a 
general rule is preserved that mode sequences are always 
from slower to faster observed drift and are concluded by 
the N mode or nulls. Nevertheless, the C mode, while not 
violating the rule, tends to occur in long PSs with little 
or no intrusion from other modes, and no examples of A 
^^ B — >■ C were found. 

In this picture the N mode (together with its accom- 
panying nulls) acts as a kind of reset, during which the 
acceleration of the driftrate ends and the subpulse pat- 
tern prepares to recommence at a slower rate, be it A, 
B or sometimes C. It is accompanied by a fading of the 
leading component and a slight shift of the central inten- 
sity to later longitudes. Interestingly, we found one "re- 
set" sequence where the leading component did remain 
strong and gave the emission a periodicity of about 12 
Pi punctuated by nulls at this periodicity. This sequence 
looked more like a fourth slow drift mode than simply a 
disordered N mode and was followed by the A mode, thus 
being consistent with the rule of incrementally increasing 
driftrates. Interestingly, the role of nulls as agents of "re- 
set" is not new: Bhattacharyya et al. (2010) noted this 
effect in B0818-41, a pulsar with variable drift, though 
not with multiple modes. 

HW comment that the modal transitions may be 
more gradual than those found in other pulsars (such 
as B0031-07). We would suggest that this may not be 
so. Initial visual inspection of the drift bands seems to 
indicate that some A — >■ B transitions may be discrete 
rather than continuous, while others may actually "over- 
shoot" the B mode, with the drift rate starting in the 
slow A mode, then quickly accelerating to even become 



8 THE POSSIBILITY OF ALIASING 

In the previous section we have analysed only the ob- 
served drift rate changes and, in common with most ob- 
servers referenced earlier, perhaps tacitly assumed that 
the observed rates are the true rate. In this section we 
analyse the consequences of challenging this assumption. 

When observing pulsars we sample the emission once 
every rotation period. If this sampling rate happens to 
be frequent compared with the timescales of intrinsic 
changes in the emission then we can obtain an adequate 
description of the phenomenon we are seeking to under- 
stand. Assuming this to be true for B1918+19, we may in- 
terpret the observed drift bands as corresponding to indi- 
vidual subbeams gradually moving across the observer's 
line of sight at the observed drift rate D (the amount of 
longitude phase the subbeam drifts in one pulse). The 
repetition rate P3 of this pattern is then P2/D and is 
significantly larger than Pi, the sampling rate. However, 
from the Irf of Fig. [S] and above all from the PS in Fig. [2] 
it is clear that P3 switches between several distinct values. 
Since P2 remains at roughly the same value, the conclu- 
sion must be that mode changes arise from changes in 
the subbeam drift rate D, which more than doubles as 
P3 reduces from mode A to mode C. 

This is a perfectly self-consistent way of interpret- 
ing our observation. However, any physical model must 
then explain why the subbeams change their intrinsic 
drift speed so drastically, and why these speeds always 
select from the same values. In the framework of the 
carousel model (Ruderman & Sutherland 1975; DR; Gil 
& Sendyk 2003) there is no suggestion as to how this can 
be achieved and why multiple fixed drift rates might be 
expected. 

Alternatively, we may assume that our sampling rate 
is relatively poor, and that we sample each subbeam at 
most twice in successive pulses as it rapidly moves from 
later to earlier phase longitude across our sightline (i.e., 
first-order aliasing). This means that the apparent drift 
rate D is no longer a measure of the true drift rate of 
each subbeam, but a confection of the much larger true 
drift D' and our sampling rate. In other words, between 
each sampling the subbeams actually drift by an amount 
comparable to the subbeam spacing P2, resulting in an 
apparent drift with P3 given by 



P3 



Pk 



(P2 - D') {Pi - 1) 



(1) 



where P3 is the true repetition rate and D' is now much 
greater than the D of the non-aliasing scenario. Now the 
P3 values measured for B1918-(-19 (6.06, 3.8 and 2.45Pi) 
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Figure 9. Here we see the same C-mode PS of Fig. [6] (bot- 
tom) folded at the assumed unaliased P!^ of 1/.59 = 1.69 Pi. 
Notice that now the drift is faster and from the traiUng to 
the leading edge, and is equivalent to a fluctuation frequency 
of -0.41 c/P\. The A- and B-modes can also be folded as in 
Fig. |6] with the appropriate values of Pg shown in Table 2, also 
giving reversed drift directions and similarity in pattern. 



are aliased quantities generated by unaliased P^ values of 
1.12, 1.36 and 1.69Pi, respectively. Again, the sequences 
of the modes can be separately folded to reveal their pat- 
terns, as in Fig. [SI but using the new P^. An example 
of this is given in Figure [51 using an unbroken sequence 
of C-mode pulses, where it can be seen the drift is from 
trailing to leading edge. 

The P3 values are well in excess of the P3 values 
which generate them because D' and P2 are close in value, 
yielding a small denominator (P2 — D') in eq. (1). Hence 
the observed, apparently large, reductions in P3 which oc- 
cur at a mode change require only a small change in D' , or 
a small change in P2, or both. Unlike the case where P3 is 
assumed to reflect the true drift, an aliased mode change 
in P3 can be observed even if the drift rate remains un- 
changed, requiring only a small, perhaps observationally 
imperceptible, change in P2. This would suggest that a 
physical account of a mode change need not be as drastic 
as under the non-aliasing assumption. It may be simply 
the result of somewhat greater or less bunching in the 
subbeams. However, the question as to why the system 
always returns to the same discrete modes remains unan- 
swered. 



9 CAROUSEL-MODEL CONFIGURATION 

The carousel model assumes that the observed subbeams 
with their corresponding P2 and D' reflect the circulation 
of a number (n) of sparks about the pulsar's magnetic 
axis. Their observed values will depend on frequency and 
height of the emission, but the values of P3, P3 and n 
should be frequency and height invariant. 



i 



In the context of B1918-fl9, the carousel model 
places an important geometric/mathematical constraint 
on P2 and D' . If the generating circle of sparks has a 
circulation time of P3 pulse periods (Pi), then the polar 
drift rate D'^^i^,. is 360° /P3 per pulse, and the separation 
of the sparks is 360° /n, where n has to be an integer. In- 
serting these into eq.(l) we obtain the intrinsic repetition 
rate to be 



P^ 



(2) 



If we now make the simplest assumption that P3 remains 
constant independently of the mode changes, then we can 
take the ratios of the mode-dependent values of P3 de- 
rived above to obtain 



riA 



0.88 



riA 



0.71 



(3) 



where ua, tlb and nc are the number of sparks (sub- 
beams) which will produce the appropriate value of P3 for 
each mode. This suggest a self-consistent model with inte- 
gral values of n where ua = 10, ub ~ 9 and nc ~ 7, and 
implies from eq.(2) a near-constant circulation time P3 
between 11.8 (from C) and 12.2 (from B), with nA = 10 
giving P3 = 12.0 exactly (a model with integral multiples 
of these values for n is theoretically valid, but unlikely to 
have physical meaning). 

What consequences does the assumption of a fixed 
circulation time of 12-Pi have on our understanding of 
the observations? From eq.(l) the observed aliased P3, 
caught by the eye and measured in the Irf, is (in units of 
A) 



Pa 



(P, - n) 



(4) 



Note that for values of n close to P3 = 12, P3 will be large 
(giving, for example, P3 = 6 for n = 10 in the A mode) 
but decrease with decreasing n. Thus we arrive at an 
interesting and self-consistent picture where the A mode 
commences with 10 subbeams, then loses a beam and 
rapidly adjusts to a B mode with 9 sparks with a slightly 
enlarged spacing, sometimes progressing to a C mode of 
7 subbeams. The spacings of the A mode beams will be 
just 70% smaller than those of the C mode. However, the 
resulting changes in P2 may be difficult to measure in 
practice (see next §) but could be a test of the model. 

Notice that this is almost precisely what can be seen 
in Fig. [4] where ffuctuations peaking near 0.83, 0.75, and 
0.58 c/Pi are evident. One can even see vestiges of a 
weak periodicity at at 0.67 c/Pi that corresponds to a 
carousel with 8 subbeams, and although we have never 
identified any stable PS with these properties, it does 
seem to be present as a transitional structure when the 
A mode sometimes "overshoots" the B and then resets 
to the B. 

The role of the N mode is interesting. We have noted 
earlier that this mode is accompanied by a periodic- 
ity with P3 « 12Pi, partly generated by repeated null 
bunches embedded in the N mode. Can it be a coinci- 
dence that this is precisely the value deduced by sepa- 
rate arguments for the fixed circulation period? There 
are two ways this might happen. Firstly, according to 
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eq.(4) this value would be consistent with a circulation of 
P3 ~ 12Pi — but with exactly n = 11 subbeams. Thus the 
sequence N -^ A ^ B ^ N corresponds to a one-by-one 
reduction in subbeams from 11 to 9 and re-set. Extrapo- 
lating further, 12 subbeams would, through aliasing, for- 
mally produce an infinite P3 — i.e., sustained non-drifting 
emission with a narrow P2 during the relatively short du- 
ration of the often irregular N mode. Thus the sequence 
may even be argued as starting at n = 12 (see Table 2). 

A second — or additional — way of explaining the N- 
mode's periodicity would be to speculate that the sub- 
beam system breaks up as n reduces below the C-mode 
value of 7. At low values of n the subbeam separation 
P2 would rapidly grow to twice or more than the values 
in the A, B and C modes (see Table 2), so that when n 
equals, say, 3 or 4 the observer would see a flickering or 
a series of pulses with deep nulls between the subbeams, 
depending on the width of the subbeams. Once the sys- 
tem has degenerated, the circulating system might just 
be partly null and partly irregularly-located beams, so 
that in the N mode we see the circulation time directly. 
Thus the N mode could be explained as both being the 
ending and the beginning of the modal cycles, empha- 
sizing its role as the mode during which the cycles are 
reset. 

Obviously missing in the sequence is the case of n = 
8 subbeams, where a P3 of 3.0 Pi would be expected. 
However, we note that there is no clear evidence of the B 
mode progressing directly to the C mode, suggesting that 
the C mode arises under separate circumstances and is 
observed to be steadier than the A and B modes (Fig. 2). 
The reason for the stability of the C mode is suggested 
by the form of eq.(4): when n is close to P3 ~ 12Pi, 
as it is in modes A and B, small modulations in P3 may 
cause observable modulations in P3 , but when n is smaller 
(7), as in mode C, modulations in the drift rate will not 
significantly modulate the aliased P3. 

A summary of this picture is given in Table 2. The 
observed values of P3, assumed to be first-order aliased, 
are shown unbracketed, together with their inferred true 
drift repetition rates P3 and inferred circulation times 
P3. The near-invariant circulation time of 12 Pi is then 
assumed to deduce the expected P3 and P3 for other in- 
tegral values of n, and these are shown in brackets. In 
the last column we give the angular spark separation for 
each n. Note how these increase weakly for high n but 
widen dramatically for small n. 

Finally, if B1918-|-19's circulation time P3 is about 
12- Pi, and the N, A, B and C modes refiect carousel 
configurations with 11, 10, 9 and 7 subbeams, respec- 
tively, then how are we to understand the 85- Pi fluctu- 
ation feature found in the Irf of in Fig. [Sf As mentioned 
above, this feature is null-associated, and while from the 
hrf of Fig. 13] it seems to be a combination of amplitude 
and phase modulation, it is carried by such small Fourier 
components that it must be a modulation of the entire 
proflle. Interesting also is the fact that one sees little of 
this feature in the analysis of shorter PSs. The A-mode 
PSs are too short, and those in the C mode do not seem 
to null. Rather it is the B and N modes that are typically 
punctuated by nulls, and the former are usually rather 



Table 2. Carousel model expectations of the circulation time 
(P3 ) for different n assuming P3 is the first alias of Pg . Com- 
parisons corresponding to observed P3 values are shown bold 
and unbracketed. 



n 


P3 


Pi 


Ps 


360/n 




{Pi) 


(Pi) 


{Pi) 


(°) 


12 


(00) 


(1.00) 


(12.00) 


30 


11 


12.0 


1.09 


12.00 


32.8 


10 


6.06 


1.20 


12.00 


36 


9 


3.8 


1.36 


12.24 


40 


8 


(3.0) 


(1.50) 


(12.00) 


42.5 


7 


2.45 


1.69 


11.83 


51.5 


6 


(2.00) 


(2.00) 


(12.00) 


60 


3 


(1.33) 


(4.00) 


(12.00) 


120 



long, averaging 53 pulses, while the N mode averages 23.5 
pulses. 

What we therefore believe, and what is evident from 
our PS in the left panel of Fig. (2] is that this feature 
is associated with the full modal cycle of N— >-A— >-B— >-N, 
which tends to start and end with nulls. Similarly, the 
43-Pi feature arises from the shorter and less frequent 
N— >-B^-N cycle. The cycle involving the C mode is of a 
different character and never exhibits this feature. 



10 EMISSION GEOMETRY 

In the previous section it has been argued on purely 
mathematical grounds that the apparent repetition pat- 
terns of B1918-|-19's emission can arise from first-order 
aliasing of equally-spaced circulating subbeams. Assum- 
ing the circulation rate to be unvarying, it was possible to 
show that each of the observed drift rates corresponded 
to a different integral number of subbeams, and hence to 
a specific angular spacing of those subbeams (Table 2). 

We must now ask whether the angular spacing gener- 
ated at the pulsar's polar cap for each mode is compatible 
with the spacing P2 of the emission beams at our observ- 
ing frequency. This spacing is only relevant for drifting 
subpulses, implying here that it must be evaluated for the 
inner cone in the central part of the profile, not the outer 
conal components on the profile edges. And, of course, we 
are interested in the values for each of the three drifting 
modes. A, B and C. 

Measuring P2 accurately is never trivial, and 
B1918-I-19 entails several further complications: first, the 
always short A-mode episodes introduce one kind of dif- 
ficulty and the incomplete trailing part of the main com- 
ponent in the C mode another. With 20/20 hindsight we 
can see that HW's P2 values in their table are inconsis- 
tent, and we will argue that they are least correct with 
their C-mode value. HW determined their values directly 
by measuring the intervals between subpulses. Another 
method uses the phase ramp of the Irf spectra, and we 
give such a diagram for the C mode in Figure [TO] There 
it can be seen both that the fiuctuation power is virtu- 
ally absent under the trailing part of the profile, but in 
the immediately adjacent leading part, the slope is such 
that a 360° excursion would correspond to about 18-20° 
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Table 3. Geometrically deduced values of n for modes A, B 
and C. 



a 


/3 


P2 


V 


n 


n 


n 


(°) 


(°) 




13.7 


-4.2 


13 


36.1 


10.0 






14.5 


40.4 


8.9 






(16) 


(44.7) 


(8.0) 






18 


50.6 


7.0 



of longitude. This gives an estimate of the C-mode P2, 
but a poor one owing to the restricted longitude range of 
the modulation. Similar phase diagrams for short B- and 
A-mode PSs (pulses 2459-2520 and the 55 A-mode pulses 
in Fig. [5| provide better estimates, suggesting P2 values 
of some 14-15° and 13±2°, corresponding to 33 and 30±5 
ms, respectively (40±4 and 33±3 ms in HW). 

P2 values can also be estimated from the hrf spectra, 
and we computed an hrf with the longitude restricted to 
the central component (not shown). This did not seem to 
significantly affect the three features corresponding to the 
drift modes. However, taking the hrf of the relatively long 
C mode PS (Figure[TT]) we can see that the primary mod- 
ulation is associated with the Fourier component peaking 
at about 20 (its harmonic very clearly peaks about com- 
ponent 40); therefore, P2 here is very near 360° /20 or 
18° or 41 ms (34±2 ms in HW). We found it difficult to 
measure the A- and B-mode P2 values by this method, 
and the problem is discernible in Fig. U The harmonic 
numbers carrying the B and A modes are clearly smaller 
as can be judged from their harmonics, but the peaks 
themselves are too broadened and conflated to give any 
reliable estimate. 

The above P2 values are given in Table 3, and we 
can use these values for each mode, together with the 
inferred geometry of the pulsar with respect to us, to 
compute what kind of carousel subbeam system underlies 
them. Starting with the profile geometry in §3 above, it 
is of importance here is that the deduced a and /3 values 
square with the central PPA rate R [= sin a/ sin /3] of - 
3.2°/°. The circumstance that this pulsar has an inside 
(poleward) sightline traverse (as does B0943-H0) is also 
quite clear from the profile geometry and linear character 
of the star's PPA traverse. 

The point is that P2 subtends a particular magnetic 
polar longitude r), which in turn can be used to estimate 
the subbeam number. DR give such a relation as their 
eq. (3), but Gupta et al.^s eq. (6) is more accurate be- 
cause it distinguishes the conal edges from the subbeam 
path within the cone. We used a version of this rela- 
tion sin 7; — sm{a + /?) sinP2/sinr where the subbeam 
path 77 is 0.8° interior to the outside conal half power 
point. Note that this value closely refiects the difference 
between the outer half-power and peak conal radii in the 
models of ET VI and Gfi & Sendyk (2000). Therefore, in 
our case, /? is negative, so the peak subbeam track has a 
magnetic polar colatitude of |/?-|-0.8°| or 3.4°. Then, the 
subbeam number n can be estimated as 360° /n. These 
simple computations are assembled in Table 3, where it 
can be confirmed that the derived subbeam numbers are 
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Figure 10. Modulation phase diagram for the C-mode PS in 
Fig.E] 
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Figure 11. The hrf for the C-mode PS in Fig. [6] Note that the 
drift modulation is carried by Fourier components around or 
just greater than 20 and that a second harmonic is discernible. 
These imply that the P2 for this mode is Ri 360° /20 = 18°. 



near integral and are fully compatible with the integers 
10, 9 and 7 for the A, B and C configurations respectively 
in Table 2. As a check, one can also carry out these com- 
putations for a positive (equatorward) sightline geometry, 
and the resultng 77 values are roughly twice as large, im- 
plying a very small number of subbeams, demonstrating 
that this sightline would be incorrect. 
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11 SUMMARY OF RESULTS 

In 1987 Hankins & Wolszczan carried out the first anal- 
ysis of the properties of B1918+19, attracted by its en- 
tirely conal triple profile and the added opportunities this 
configuration presents to investigate the relationships of 
emission in its inner and outer cones. Using Arecibo PSs 
they were able to identify the four basic single-pulse 
modes of emission (A, B, C and N) and noted specific 
cyclic relationships between these modes. 

Here, with the advantage of more recent data from 
the same telescope but at a slightly lower observing fre- 
quency, we explore the geometry of the pulsar and re- 
fine and add to the known single pulse repetition pat- 
terns. Further, we propose a possible carousel configu- 
ration present at the polar cap which may explain the 
mode changes with a simple idea. Before discussing this 
model in the next section, it is useful to summarise our 
observational results; 

• B1918-fl9 has integrated properties common to 
many pulsars: it has a conal triple (cT) profile with very 
nearly the usual angular dimensions for its inner and 
outer cones — half-power conal radii p at 1-GHz of 4.1 and 
6.5° — as computed originally in ET VI. Its integrated po- 
larisation properties suggest that the magnetic latitude 
a and sightline impact angle /3 are some 13.7 and -4.2°, 
respectively, such that its inside (poleward) sightline tra- 
verse cuts its inner cone tangentially (/3/p ~ —0.90) and 
outer cone obliquely (/3/p ~ —0.65) 

• The profile exhibits leading/trailing asymmetry of 
two types: the leading (outer) conal component is much 
stronger than the trailing one; and the leading part of 
the main (inner cone) component is stronger than the 
trailing part. This asymmetry is least marked in the N 
mode, and becomes stronger in the A and B as well as the 
C mode — wherein trailing portions of the profile nearly 
vanishes. 

• As HW found, the star's pulse sequences exhibit 
three highly discrete drifting modes. A, B and C, as well 
as a non-drifting N mode. However, the pulsar exhibits 
five prominent fiuctuation features, three corresponding 
to the A, B and C drift modes with P3 values near 6.0, 
4.0 and 2.5 Pi, and two null-related features at about 12 
and 85 Pi. 

• About 9% of the pulses are "nulls" but these nulls are 
distributed in a highly non-random manner, evidenced 
by the substantial fraction of nulls with long (2-10 pulse) 
durations. The A and C modes seem not to null, and the 
B mode infrequently, so the nulls occur mainly within the 
N mode at a roughly 43% fraction. 

• Fluctuations in the inner-cone main and outer-cone 
leading components are phase-locked in all three modes — 
and with very similar phase relationships. This suggests 
that the drift features. A, B and C, may be seen as con- 
sisting of a similar pattern across both inner and outer 
cones spun past our sightline at three different speeds, or 
as a widening pattern at a common speed. The phase- 
locking is very reminiscent of the double-cone model pro- 
posed for B0818-41 by Bhattacharyya et al. (2007), and 
has been argued for B1237+25 (Srostlik & Rankin 2005), 
B0836-24 (Gupta et al. 2004, Esamdin et al. 2005) and 
B1039-19 (Bhattacharyya et al. 2011). Finding this same 



property in B1918-f 19 mode-by-mode suggests that lock- 
ing between inner and outer cones is a common, maybe 
ubiquitous, feature in double-cone pulsars with regular 
drift. 



12 CONCLUDING REMARKS 

For B1918+19 the classic carousel model of Ruderman 
& Sutherland (1975) predicts a circulation time P3 of 8.5 
Pi, well below what would be expected if the observed P3 
values were interpreted directly as measures of spark-to- 
spark transition speeds. With a more subtle model involv- 
ing a central discharge (Gil & Sendyk 2000) it may well 
be possible to achieve a more plausible circulation speed, 
but it remains hard to understand why the circulation 
time switches between three fixed but very different val- 
ues in this and other pulsars (see, for example, discussion 
of the quantised drift patterns of B2319-(-60 in Wright & 
Fowler 1981). 

Here we have explored the view that the observed 
driftbands result from first-order aliasing of a much faster 
drift. This is what the original carousel model (assuming 
any reasonable number of "sparks") would actually pre- 
dict, since our sampling rate (Pi) is comparable to the 
frequency (P3) with which the "sparks" are presented to 
u^j. This is not the first time subpulse drift has been as- 
cribed to aliasing {e.g., DR, Janssen & van Leeuwen 2004, 
Gupta et al. 2004, Bhattacharya et al. 2007), but here we 
extend the idea to account for the presence of quantised 
drift rates and their sequence. Using our observations we 
are able to deduce that aliasing requires 10, 9 and 7 sub- 
beams corresponding to the A, B and C modes and, more 
surprisingly, that these are all consistent with a common 
circulation time of about 12 Pi, not far from the orig- 
inal RS figure. This essentially mathematical-geometric 
picture appears to be compatible with the observed P2 
values, which increase gradually as P3 falls, as the model 
requires (see Table 2). Thus we no longer need to ex- 
plain a changing circulation time, and can understand 
the changing — quantised — drift rates as a progressive re- 
duction in the number of sparks. This is physically much 
easier to understand. 

It is significant that the N mode is also clearly linked 
to the circulation time. The N mode was not involved in 
the deduction of A, B and C's common P3, and yet its 
repetition rate P3 ~ 12 - discovered here - is virtually 
the same as P3. So what does this mode represent? We 
demonstrate in §9 that the mode can betray the under- 
lying circulation time by firstly ending the mode cycle 
through the vanishing subbeams with wide angular sep- 
arations, and then by re-starting the mode cycle through 
the aliasing of 11 subbeams. Thus the mode sees the re- 
generation of the mode cycle and must be the phase in 
that cycle when the subbeams re-establish themselves. 

One may speculate whether carousel configurations 

•^ Note that higher orders of aliasing are perfectly possible in 
mathematical terms. These will require even faster circulation 
times and have not been preferred here for physical reasons. 
However, they cannot be excluded for different physical con- 
figurations and in carousels of other pulsars. 
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with varying numbers of subbeams might explain the 
multiple different drift modes exhibited by other pul- 
sars. The best-known examples, B0031--07 (Huguenin et 
al. 1970; Vivekanand & Joshi 1997), B2319+60 (Wright 
& Fowler 1981), B2303+30 (Redman et al. 2005) and 
B1944-I-17 (Deich et al. 1986, Kloumann & Rankin 
2010), have similar ratios in their Pi values to those of 
B1918-I-19 and are worthy of further investigation in the 
context of aliased drift. 

However, as is evident in the sequences of Fig. 2 and 
as noted above, the drift rates of B1918-I-19 are not to- 
tally steady in the A and B modes and at the transition 
can initially overshoot their new value. This suggests that 
the circulation time is not held rigidly steady as supposed 
in the basic calculations here. Pulsars B2016-I-28, B0818- 
41 and B0826-34 all intermittently display curved drift- 
bands, so changes in pulsar drift patterns cannot always 
be explained by the sudden disappearance of a subbeam. 

The point we make in this paper is that subtle vari- 
ations in the circulation time and in the number of circu- 
lating beams can, if the sampling rate of once per pulse 
period is infrequent for probing the true circulation time, 
generate the extraordinarily dramatic and spectacular ef- 
fects seen in the pulse sequences of many pulsars. 
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